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Abstract 

We perform a lattice QCD study of the p meson decay from the Nf = 2 + 1 full QCD config- 
urations generated with a renormalization group improved gauge action and a non-perturbatively 
0(a)-improved Wilson fermion action. The resonance parameters, the effective p — > tttt coupling 
constant and the resonance mass, are estimated from the P-wave scattering phase shift for the 
isospin 1 = 1 two-pion system. The finite size formulas are employed to calculate the phase shift 
from the energy on the lattice. Our calculations are carried out at two quark masses, m n = 410 MeV 
(m^/mp = 0.46) and m n = 300 MeV (m^/mp = 0.35), on a 32 3 x 64 (La = 2.9 fm) lattice at the 
lattice spacing a = 0.091 fm. We compare our results at these two quark masses with those given 
in the previous works using Nf = 2 full QCD configurations and the experiment. 

PACS numbers: 12.38. Gc, 11.15.Ha 
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I. INTRODUCTION 



Recent progress of simulation algorithms, supported by the development of computer 
power, has made it possible to study hadron physics at the physical quark mass by lattice 
QCD (see Ref. [l| for recent reviews), and lattice calculations have clarified the properties of 
many hadrons. The studies are mostly concentrated on stable hadrons, however. Resonances 
pose an important issue both in terms of methodologies and physical results. 

Among the resonances, the p meson is an ideal case for the lattice calculations, because 
the final state of the decay is the two-pion state which can be treated on the lattice pre- 
cisely. In the early stage of studies of the p meson decay, the transition amplitude (vr7r|p) 
extracted from the time behavior of the correlation function (7i(t)7i(t)p(0)) was used to 
estimate the decay width, assuming that the hadron interaction is small enough so that 
(mr\p) < | (p|p> (ttttItttt) I 1 / 2 is satisfied 043]- 

A more realistic approach is a study from the P-wave scattering phase shift for the 
isospin 1 = 1 two-pion system. The finite size formulas presented by Liischer in the center 
of mass frame 0] and extensions to non-zero total momentum frames 0, 0] are employed 
for an estimation of the phase shift from an eigenvalue of the energy on the lattice. The 
first study of this approach was carried out by CP-PACS Collaboration using Nf = 2 full 
QCD configurations (m n = 330 MeV, a = 0.21 fm, La = 2.5 fm) [9]. After this work ETMC 
Collaboration presented results with Nf = 2 configurations at several quark masses (m 7 
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290, 330 MeV (La = 2.5 fm), m w = 420, 480 MeV (La = 1.9 fm), a = 0.079 fm) 
Recently Lang et al. reported results of high statistical calculations with the Laplacian 
Heaviside smearing operators on a single Nf = 2 gauge ensemble (m n = 266 MeV, a = 
0.124 fm, La = 1.98 fm) [l2 . 

In the present work we extend these studies by employing Nf = 2 + 1 full QCD configura- 
tions and working on a larger lattice volume. Our calculations are carried out with the gauge 
configurations previously generated by PACS-CS Collaboration with a renormalization group 
improved gauge action and a non-perturbatively 0(a)-improved Wilson fermion action at 
/3 = 1.9 on 32 3 x 64 lattice (a = 0.091 fm, La = 2.9 fm) [13]. We choose two subsets of the 
PACS-CS configurations. One of them corresponds to the hopping parameters K u d = 0.13754 
for the degenerate up and down quarks and k s = 0.13640 for the strange quark, for which 
the pion mass takes fn n — 410 MeV [vn^ j viip — 0.46). The other is at n u d — 0.13770 and 
k s = 0.13640, corresponding to = 300 MeV (m^/rrip = 0.35). We extract the scattering 
phase shift on three momentum frames, the center of mass and the non-zero momentum 
frames with the total momentum P = (27r/L)(0, 0, 1) and P = (2n/L)(l, 1, 0), to obtain the 
phase shifts at various energies from a single full QCD ensemble as in the previous works 
by ETMC [Io|, G3| and Lang et al. Q. 

We note that QCDSF Collaboration calculated the scattering phase shifts for the ground 
state in the center of mass frame at several quark masses, (m^ = 240 — 430 MeV) [13]. They 
estimated the resonance parameters from these results, assuming that the effective p — > titt 
coupling constant does not depend on the quark mass. BMW Collaboration presented their 
first preliminary results with Nf = 2+1 configurations (m, = 200, 340 MeV, a = 0.116 fm) at 
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Lattice 2010 15j. We also refer to works exploring an application of the stochastic Laplacian 
Heaviside smearing to the two-pion states with the isospin / = 0,1,2 using Nf = 2 + 1 
configurations on the large lattice volume in Ref. [l6j |. 

This paper is organized as follows. In Sec. |TT]we give the method of the calculation. The 
simulation parameters of the present work are also presented. We present our results and 
compare ours with those by the other works in Sec. IIII1 Conclusions of the present work are 
given in Sec. IIV1 Result of a pilot study of the present work at m n = 410 MeV was reported 
at Lattice 2010 (l7j . All calculations are carried out on the PACS-CS computer at Center 
for Computational Sciences, University of Tsukuba. 



II. METHODS 

A. Finite size formula 

In order to calculate the P-wave scattering phase shift for the isospin 1=1 two-pion 
system at various energies from a single full QCD ensemble, we consider three momentum 
frames, the center of mass frame (CMF), the non-zero momentum frames with total mo- 
mentum P = (2tt/L)(0,0, 1) (MF1) and P = LJ 3) (MF2). These frames have 
been also considered in the previous works by ETMC [lfj, lll| and Lang et al. 12| . In these 
momentum frames the P-wave state is decomposed as 



(1) 



where P is the total momentum, g is the rotational group in each momentum frame on the 
lattice and V is the irreducible representation of the rotational group. In the present work 
we consider four irreducible representations : Tj~ in the CMF, A^~ and E~ in the MF1, and 
in the MF2. The ground and the first excited states of these representations with the 
isospin (I,I Z ) = (1,0), ignoring the hadron interactions, are shown in Table [H 

The scattering phase shift is related to an eigenvalue of the energy on the lattice by 
the finite size formula. The formula in the CMF was presented by Liischer jsj , that in the 
MF1 by Rummukainen and Gottlieb Q, and the MF2 by ETMC Q. The formulas for the 
representations considered in the present work are written by 
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(2) 



V 00 for in CMF 

V 00 - V 20 for E~ in MF1 

tan5(fc) ~ ) V 00 + 2 V 20 ^ for A^ in MF1 ' 

V 00 - V 20 + VQV 22 for in MF2 

with the P-wave scattering phase shift 5(k). The real function V/ m is defined by 

V lm (k; P) = j- , e -*W4 . z ( V q d ) ( 3 ) 



3 



with d = PL/(2ir) and q = kL/(2ir), where P is the total momentum and k is the two- 
pion scattering momentum defined from the invariant mass \/s by y/s = \J E 2 — |P| 2 = 
2 \Jk 2 + m 2 with the energy E in the non-zero momentum frame. In ([3]) 7 is the Lorentz 
boost factor from the non-zero momentum frame to the center of mass frame given by 
7 = E/sfs. The function Zi m (l; q; d) in ([3]) is an analytic continuation of 

Z lm (s;q;d)= ^ m (r)-(|r| 2 -g 2 )- s , (4) 

reD(d) 

which is defined for Re(s) > (7 + 3)/2, where 3^m( r ) is a polynomial related to the spher- 
ical harmonics through 3^m(r) = ■ Yi m (Q) with Q the spherical coordinate for r. The 
convention of Yi m (Q) is that of [18|. The summation for r in (J4]) runs over the set, 

D(d) = { r I r = ^(n + d/2) , n G Z 3 } . (5) 

The operation 7 _1 is the inverse Lorentz transformation : 7 _1 x = X11/7 + x_|_ , where 
X|| = (x-d)d/d 2 is the parallel component and xj_ = x — X|| is the perpendicular component 
of the vector x in the direction d. Zi m (l; q; d) can be evaluated by the method described in 



Ref. 191. 



B. Extraction of energies 

In Fig. [1] we show values of the invariant mass \fs divided by the p meson mass m p 
for the states tabulated in Table fl] on our gauge configurations at m n = 410 MeV (upper 
panel) and m n = 300 MeV (lower panel). Here we ignore the hadron interactions and use 



the values of m p and m n obtained in the previous work in Ref. [13J . For the Tj" and the E~ 
representation, we only calculate the scattering phase for the ground state in the present 
work. From Fig. [1] the energies of these states are expected to be much smaller than those 
of the excited states, even if the hadron interaction is switched on. Thus, we extract the 
energies of these states by a single exponential fit for the time correlation functions of the 
p meson as carried out in a usual study of the hadron spectrum. We use the local p meson 
operator for the sink and a smeared operator for the source as discussed later. 

For the A^" and the Bj" representation, we also calculate the scattering phase shift for 
the first excited state. In order to extract the energies of the lower two states for these 



representations, we use the variational method 20j with a matrix of the time correlation 
function, 

rm = ( <0| (7T7r)t(t)M(g |0> (0| (7nr)t(t) p(ts) |0> \ (6] 
U ^(01^(^(^)^)10) (0\p^t)p(t s )\0) J ' U 

for each representation. The energies are extracted from two eigenvalues X n (t) (n — 1, 2) of 
the matrix, 

M(t) = G(t)G-\t R ) , (7) 

with some reference time £r. 
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Here we comment on the discussion of the generalized eigenvalue problem (GEVP) by 
ALPHA Collaboration in Ref. 21] . In a 2 x 2 matrix case of G(t), they proved that the 
effective mass of the eigenvalue X n (t) (n = 1, 2) of the matrix M(t) in (J7]) can be written as 



E n + 0{e- (Em ~ En)t ) (n = 1, 2 , m > 3) , (8) 

in a large time region with the eigenvalue of the energy Ei (i = 1, 2, • • • ). Here it should be 
noted that their proof was given only for the case where G(t) is a Hermitian matrix. In our 
case we use different operators for the sink and source in G(t) as explained later. Therefore 
G(t) is not a Hermitian matrix and the discussion of GEVP cannot be applied to our case. 
In the present work we assume that the lower two states dominate G(t) in a large time 
region. This is expected to be a good approximation, because the second excited state takes 
a much higher value ( E 3 /E 2 > 1.3 in the absence of hadron interactions for both quark 
masses studied in the present work ). With this assumption, the second term of (jSJ) does 
not appear for a general case of G(t), and the energy for the ground and the first excited 
states can be extracted by a single exponential fit for the eigenvalue X n (t) (n = 1,2) in a 
large t region. 

In (jBJ) the operator p(t) is given by 

3 

p(t) = J2PrPj(P^)/\p\ » (9) 

where Pj(p,t) is the local operator for the neutral p meson at the time slice t with the 
momentum p. The momentum takes p = (2tt/L)(0, 0, 1) for the A 2 and p = (2ir/L)(l, 1, 0) 
for the B]~ representation. Hereafter we assume that the momentum p takes one of these 
two values depending on the representation. 

In ([6]) (iT7r)(t) is an operator for the two pions with the momentum and p, which is 
defined by 

(7T7r)(t) = i=(7r + (0,t 1 )7r-( P ,t)-7r-(0,t 1 )7r + ( P ,t)) xe" 1 -^"*) , (10) 

where vr ± (p, t) is the local pion operator with the momentum p at the time slice t. The time 
slice of the pion with the zero momentum is fixed at t\ 3> t, and the time slice of the other 
pion t runs over the whole time extent. An exponential time factor in (fTUj) is introduced so 
that the operator has the same time behavior as that of the usual Heisenberg operator, i.e., 

(0| (tttt)^) = (0| (7r7r)t(0) e~ m for t x > t , (11) 

with the Hamiltonian H. 

In the previous works the time slices of the two operators for the pion in the sink operator 
( flOj) are set equal t\ = t, and they simultaneously run over some time interval. In that case 
we need to repeat solving quark propagators for the time slices in that time interval. This 
computer-time consuming procedure can be avoided by fixing the time slice of one of the 



5 



pion at t\ as done in the present work. But we need to set t\ 3> t to avoid contamination 
from higher energy states produced by the operator at ti in this method. 

Two operators (7T7t)(£ s ) and p(t s ) are used for the sources in ([6]), which are given by 

(W)(t s ) = (7T + (0, t s )7T-( P , t s ) - 7T-(0, t s )7T + ( P , t a )) , (12) 

P^') = l^S^f^^p^t.) -S( Z ,t s )7^(z,t s ))e ip ' z , (13) 

where 7 P = Y^=iPj ' Tj'/IpI- The operator Q(z, t s ) (Q = U, D) is a smeared operator for the 
up or the down quark given by 

Q(z,t,) = £; g (x,t,)-*(|x-z|) , (14) 

x 

where g(x, t s ) (q = u, d) is the up or the down quark operator at the position x and the time 
t s . We adopt the same smearing function \l/(|x|) as in Ref. [13J, i.e., ^(|x|) = Aexp(— B\x\) 
with \I/(0) = 1 and the smearing parameters : (A,B) = (1.2,0.17) for m n = 410 MeV 
and (A, B) = (1.2,0.09) for m n = 300 MeV. The operator (|T4j) is used after fixing gauge 
configurations to the Coulomb gauge, assuming that an ambiguity of the gauge fixing does 
not cause a significant systematic error in the time correlation function. In ([TBI a summation 
over z is taken to reduce a statistical error and 

T = { z | z = (L/2) • (m, n 2 , n 3 ) , rij = or 1 , iV r = 8 } (15) 

is chosen in the present work. The smeared operator ffTB"]) is also used to extract the energy 
of the ground state for the Tj" and the E~ representation, setting the momentum p = 
and p = (27r/L)(0, 0, 1), respectively. 

Here we note that the operator (inr) (t) in ([TO]) is not an eigenstate under exchange of the 
momenta of the two pions. Thus it has no definite parity and does not belong to the irre- 
ducible representation of the for p = (27r/L)(0, 0, 1) and the D 2 h for p = (2tt/L)(1, 1, 0). 
It is actually a linear combination of components of two irreducible representations A2 + 
for the D^h arid + A] 1- for the -D2/1. The other operators in G(t) in ([6]), however, belong 
to the irreducible representation Ag or Bj~ depending on the momentum p. Therefore, it 
is not necessary to worry about mixing to other irreducible representations in G(t) in ([6]). 



C. Calculation of G(t) 



The quark contractions of the components of the matrix of the time correlation function 
G(t) in (EJ) are shown in Fig. |2j The time runs upward in the diagrams. The vertices refer 
to the pion or the p meson operator with the momentum at the time slice specified in the 
diagrams. The momentum p takes p = (27r/L)(0, 0, 1) for the A^" and p = (2ir/L)(l, 1,0) 
for the B]" representation. The p meson operators at t s are the smeared operators and the 
other is the local operator. 
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We calculate the quark contractions in Fig. [2] by the source method and the stochastic 
noise method as in the previous work by CP-PACS [9]. We introduce a U(l) noise £j(x) 
which satisfies 



(16) 



£j(x)£;(y) = <5 3 (x - y) for N R -> oo , 

3=1 

where iV^ is the number of noises. We calculate the following four types of quark propagators: 



Q AB (x,t\q,t s ,Q = ^(D-^ABfatwt,) ■ [e iq %(y 



(17) 



W AB (x, t\k, t a \q, t s , Q = Y] V( J D- 1 ) AC (x, t; y, t a ) ■ [e ik ^ 75 Q(y, t a |q, t fl , (18) 



y c 



Q AB (x, t\z, t s ) = J2( D ~ 1 )ab(*, t- y, t s ) ■ [tf (|y - z| 
y 

W AB {x,t\k,t a \z,t s ) =^2^2(D- 1 ) AC (x,t;y,t a ) 

y c 



e lk - y 7 5 Q(y,ta|z,t s ) 



CB 



(19) 
(20) 



where A, B and C refer to color and spin indices, and ^(|x|) is the smearing function in 
(fl4|) . The square bracket is used as the source for the inversion of the Dirac operator D. 
The propagators for the smeared quarks (Q and W) are solved on the gauge configurations 
fixed to the Coulomb gauge, while the gauge is not fixed for those of the stochastic quarks 
(Q and W). 

The function G nn ^. nw (t) for the first diagram in Fig. [2] can be calculated by introducing 
an another U(l) noise %(x) having the same property as £j(x) in (Fl6|) . 

N R 

GSSnnrCt) = ^^(g t (x,t 1 |0,t s ,^)g(x,t 1 |0,t s ,r /j ))e m -(' 1 ^ 

3=1 x,y 

x e-^ y (gt(y ) i|o ) t s ^ i )g(y ) t| P)is) ^.)^ , (21) 

where the bracket means trace for the color and the spin indices. The exponential time 
factor comes from the definition of the operator for the two pions in (fit)]) . The G nn ^. nn (t) 
for the second diagram is given by exchanging the momentum and the time slice of the sink 
in d2H). 

The function G^^^t) for the 3rd to 6th diagrams can be obtained by 

Nr 

GfXAt) = ^^e-^ x (^t (X)t |o ) t 1 |o,t s ,^0^(x,t|O,t,|p,t s ,e i ))e m -( t1 ^ , 

3=1 x 

Nr 

GjSLr(t) = ^^e-p- x (^t (X)t |o, tl |o,t s ,e J o^(x,t|p,t s |o,t s ,e i ))e m - {ti - t) , 



3=1 x 

Nr 



GfX,(t) = ^^e-^^^tlCtxIC^O) W\x,t\0,t s \ -p,t. l ^))e m "' ( * 1 ^ 

3=1 x 
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N R 

GfXM = ^^e- ip - x (^(x,t|0,t 1 |0,t s ,0)^ t (x,t|-p,t,|0,t s ,0))e m - (tl - t) (22) 

j = l X 

The function G n7T ^ p (t) for two diagrams in Fig. [2] can be similarly calculated by 

N R 

Ggft p (f) = ^^e- l P- x (^t (X)t |_ P)ts | 0) t S) ^)(757 P )Q(x,t|0,^0)) , 

j = l X 

G [ £%(t) = ^^e- l P- x (gt( x ,t|o,t s ,e,) ( 75 7 P ) W(x,t\p,t s \0,t s ,Q) , (23) 

j = l X 

where 7 P = Y%=iVj • Tj/IpI- 

We can obtain the function Gp^^t) for two diagrams in Fig. |2]by 

r zer x 

G^i(t) = X) X) e-^-^^W^Cx, t|0, *.) *.) (w P ))e m ^-» (24) 

r zer x 

where T and N? are defined by (|T5l) . The component G p ^ p (t) is given by the Q-type 
propagators as the usual time correlation function for the p meson, 

<W*) = -^^e i P- z ^e-^ x (g t (x,t|z,t s )( 75 7 P )Q(x,t|z ! t s ) ( 7s 7 P )) • (25) 
r zer x 

We calculate the Q-type propagators (IT7|) for combinations of q and U(i) noise : 

(q, noise) = { (0, £), (0, r/), (p, 0, (~P, } > (26) 

and VT-type propagators ( 1T81) for combinations of k, i a and q : 

(k, t a \q) = { (p, t s \0), (-P, * s |0), (0, f s |p), (0, t s \ - p), (0, tx|0) } , (27) 

using the same U(l) noise £ in common. The Q-type ( 1T91) and the W-type propagator ()20l) 
for (k, t a ) = (0, ii) are solved for the set z G T. 

D. Simulation parameters 

Calculations in the present work employ Nf = 2 + 1 full QCD configurations previously 
generated by PACS-CS using a renormalization group improved gauge action and a non- 
perturbatively 0(a)-improved Wilson fermion action at /3 = 1.9 on 32 3 x 64 lattice (a = 
0.091 fm, La = 2.9 fm) [13]. We choose two subsets of the PACS-CS configurations. One of 
them corresponds to the hopping parameters K ud = 0.13754 for the degenerate up and down 
quarks and k s = 0.13640 for the strange quark, for which the pion mass takes m n = 410 MeV 
{m^/rrip = 0.46). The total number of configurations analyzed every 10 trajectories is 440. 
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We estimate the statistical errors by the jackknife method with bins of 400 trajectories. 
The other set is at K u d = 0.13770 and k s = 0.13640, corresponding to m n = 300 MeV 
(m 7T /m p = 0.35). The total number of configurations of this set is 400 and the measurements 
are done every 20 trajectories. The statistical errors are estimated by the jackknife method 
with bins of 800 trajectories. 

The periodic boundary conditions are imposed for both spatial and temporal directions 
in configuration generations. We impose the Dirichlet boundary condition for the temporal 
direction at t — and t = T in calculations of the quark propagators, to avoid the unwanted 
thermal contributions produced by propagating two pions in opposite directions in a time. 
For both quark masses, we set the source operators (im)(t s ) in (fT2|) and p(t s ) in ( TT3]) at 
t s = 12 to avoid effects from the temporal boundary, and the zero momentum pion in the 
sink operator (inr)(t) in ( ITUj) at ti = 42. 

In order to see effects from the Dirichlet boundary, we calculate the time correlation func- 
tions of the single meson channels setting the source at t s = 12 with the Dirichlet boundary 
condition. We compare these with those obtained with the periodic boundary condition in 
the previous work in Ref. 13|] , and find no difference beyond statistical fluctuations between 
them. This shows that t s = 12 is sufficiently large to avoid effects from the boundary for 
the single meson channels. We assume from this that t s = 12 and t\ = 42 (22 away from the 
boundary) are safe distances to avoid effects from the boundary also for the time correla- 
tion functions of the two-pion state. We find that the effective mass of the time correlation 
function for the pion with the zero momentum reaches plateau after t — t s = 12 for the both 
quark masses. We can expect from this that the eigenvalues X n (t) (n = 1,2) of the matrix 
M(t) in ([7]) have a single exponential behaviors in a time range t < t\ — 12 (= 30). 

In a pilot study of the present work at = 410 MeV, calculations of the phase shift 
were carried out only for three representations, Tj~, E~ and , setting the number of 
the random noise Nr = 10 in (fl6l) . The results on this pilot study have been reported 
in Ref. (l7| . We found from this study that errors from a finiteness of the number of the 
random noise are small enough compared with the statistical error for three representations 
even if we take Nr = 2. Assuming that this observation applies also for both quark masses 
and all representations in the present work, we set Nr = 2 for all calculations after this pilot 
study. 

In order to reduce the statistical error, we carry out additional measurements shifting 
the time slice of the source operators t s , the zero momentum pion in the sink operator 
t\ and the Dirichlet boundary condition by a time shift At simultaneously, and average 
over these measurements. For = 410 MeV, the measurement of the pilot study for the 
three representations was done without shifting the time slices. We carry out additional 
measurements for all representations with the time shift At = T/2 and T/4. We average 
over these two measurements for the B]~ representation and all three measurements for 
three representations, Tj~, E~ and Ag. For m n = 300 MeV, the measurements for all 
representations are carried out with the time shift At = and T/2, and are averaged. 
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III. RESULTS 



A. Time correlation function 

In Fig. E] we show the real part of the diagonal components (tttt — > tttt and p — > p ) and 
imaginary part of the off-diagonal components (p — > tttt and tttt — > p) of the matrix of the 
time correlation function G(t) in ([6]) for the and the Bj~ representation at m n = 410 MeV. 
We note that the diagonal components are real and the off-diagonal components are pure 
imaginary by P and CP symmetry. Choosing = 23 as the reference time of the variational 
method for the matrix M{t) in ([7]), we obtain the two eigenvalues Ai(t) and \2{t) of the 
matrix, which corresponds to the correlation function for the ground and the first excited 
state respectively, for each representation. 

The effective masses of the time correlation functions for six states considered in the 
present work at m n = 410 MeV are plotted in Fig. HJ We can find plateaus in the time 
region t > 23. The results of the energies E extracted by a single exponential fit for these 
time correlation functions are tabulated in Table El together with adopted fitting ranges. 
We choose smaller value for the maximum time of the fitting range for the and the 
representation than those for the others to avoid contamination from higher energy states 
produced by the zero momentum pion at t\ = 42 in the sink operator (7T7r)(t) in (TlPl) . In 
Fig. HI the results of the fitting with one standard deviation error band are also expressed 
by solid lines. The dotted line for the A^~ and B]~ representation in the figure indicates the 
energy of the two free pions for each representation. 

The components of the matrix of the time correlation function G(t) at m n = 300 MeV 
are plotted in Fig. [5] and the effective masses in Fig. [BJ where tn = 23 is also chosen as the 
reference time. The statistics is less than that at = 410 MeV, but we also see plateaus 
in the effective masses for t > 23. The results of the energies E extracted by a single 
exponential fit are tabulated in Table IIHI 

In the previous work by CP-PACS, carried out at the lattice spacing a = 0.21 fm, they 
found a large violation of the continuum dispersion relation for the single pion state due to 
the discretization error on their gauge configurations. The discretization error also affects 
calculations of the invariant mass a/s and the scattering momentum k for the two-pion 
system since they are evaluated from the energy E. The continuum relation is given by 
\/i = \J E 2 — |P| 2 = 2 a/ k 2 + m 2 , while there are several alternatives on the lattice, e.g., 



The two-pion scattering momentum cannot be uniquely defined due to the breaking of the 
translational and rotational symmetries in the finite lattice spacing as mentioned in Ref. 
The momentum k given by (f29j) is just one of the choices of the momentum, thus the 
discretization error cannot be fully avoided by using ( )28l) and ( 1291) . In the work by CP- 



3 




(28) 



2 ■ sin 2 (A;/2) = cosh(v/s/2) - cosh(m^) . 



(29) 
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PACS, they regarded the difference of the final results for the choice of the relations as the 
systematic error from the discretization error. 

We also monitor the validity of the continuum dispersion relation for the single pion 
state and find that the violation is negligible in the present work, one reason being that 
our lattice spacing a = 0.091 fm is much smaller than that for the CP-PACS case. We 
compare the energy E extracted from the the time correlation function with that given by 
the dispersion relation E e g = a/|p| 2 + mf from the mass and the momentum p. The 
results for = 410 MeV are E/E eS = 0.9988(15) for p = (2tt/L)(0, 0, 1) and 0.9988(63) 
for p = (2tt/L)(1, 1,0). Those for m n = 300 MeV are 0.9924(81) and 0.984(16). Therefore, 
the violation of the dispersion relation for the single pion state is not seen on our gauge 
configurations. We also evaluate y/s for the two-pion state by ( 128]) , but we see no difference 
over the statistics from that given by the relation in the continuum. From this we calculate 
\fs and k from E by the relation in the continuum, avoiding ambiguities possibly caused 
by the choice of the relations. The results given by this way are tabulated in Table Hi] and 
Table ED 



B. Scattering phase shift and resonance parameters 

The scattering phase shift 5(k) obtained by substituting the scattering momentum k and 
the total momentum P into the finite size formulas in ([2]) are presented in the lower part 
of Table [TT1 for m n = 410 MeV and Table [ED for m w = 300 MeV. We use the lattice spacing 



determined from rriQ in Ref. [13[, a = 0.0907(13) fm (1/a = 2.176(31) GeV), to get the values 
in the physical unit, where the error of the lattice spacing is not included. In Fig. [7] the 
results of fc 3 /(-y/s • tan5(&)) are plotted as a function of square of the invariant mass s for 
= 410 MeV (upper panel) and m n = 300 MeV (lower panel). The finite size formulas 
for the A 2 and the Bj~ representation are plotted by dotted lines. Divergent point on these 
lines corresponds to the square of the invariant mass of the two free pions. In the figure 
the error bars of s and fc 3 /(-y/s • tan5(/c)) are plotted regarding them as independent. But 
these are fully correlated by the finite size formula, so the true error lies along the dotted 
line corresponding to the formula. Thus, the error bars in the figure indicate projections of 
the true error bar on the finite size formula to the vertical and the horizontal axis. 

In order to extract the resonance parameters from the results of the scattering phase 
shift, we try to parametrize the resonant behavior of the P-wave phase shift in terms of the 
effective p — > nn coupling constant g p7T7T as 

l " - 6r K - s) , (30) 



tan5(fc) g*^ p 

where m p is the resonance mass and g p7Tn is defined though the effective p — > nn Lagrangian 

as 

L c k = g pTYir ^2 e ^c(h - k 2 ) tl p a ^(p)7f b (k 1 )TT c (k 2 ) . (31) 

liabc 
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This parametrization has been widely used in the previous works of the p meson decay. The 
p meson decay width at the physical quark mass is related to the coupling constant by 

Tp = ^^W- = 4.237 MeV x , (32) 

where mJJ h = 775.5 MeV is the actual p meson mass and (fc ph ) 2 = (mP h ) 2 /4 — (m^ 11 ) 2 
(mf = 135 MeV). 

By chi-square fitting of the scattering phase shifts with the fit function ( |30l) . we obtain, 



gfmv = 5.52 ± 0.40 , 
am p = 0.4103 ± 0.0026 , 
m p = 892.8 ± 5.5 ± 13 MeV , (33) 



for m n = 410 MeV, where the first error of m p is the statistical and the second is the 
systematic uncertainty for the determination of the lattice spacing. In the fitting, we define 
the chi-square for each data point by squaring the ratio of the distance from the data point to 
the fitting line (I30j) along the finite size formula and the true statistical error calculated along 
the finite size formula. The errors of the resonance parameters g pn7T and m p are estimated 
by the jackknife method as for the other values. In the upper panel of Fig. [7] we draw a 
fitting line by a solid line. We can find that the fitting with the function ( 1301 goes well in 
the large energy region at = 410 MeV. 

For = 300 MeV the statistics of our data is not enough to discuss a quality of the 
fitting with the fit function ([30]) as shown in Fig. [7J Improving the statistic by using some 
efficient smearing techniques for the two-pion operator may be necessary for an investigation 
of a reliability of ( 1301) . We must leave this issue to studies in the future. Here we carry out 
the chi-square fitting as done at = 410 MeV, assuming that the function (|30|) also works 
well in our energy region at m n = 300 MeV. The results of the fitting are given by 

g p7T7T = 5.98 ± 0.56 , 
am p = 0.396 ±0.010 , 
m p = 863 ± 23 ± 12 MeV , (34) 

where the second error of m p is the systematic uncertainty for the determination of the 
lattice spacing. We draw a fitting line by a solid line in the lower panel of Fig. [71 

From (|33D and ([34]) we find that the g p7T7T at the two quark masses are consistent within 
the statistical error and also with the experiment g p7Tn = 5.874 ± 0.014 given from the 



experimental result of the decay width T p = 146.2 ± 0.7 MeV [22] by ( |32l) . This suggests 
a weak quark mass dependence of the coupling constant. But our calculations are carried 
out only at the two quark masses and a reliability of ( 130]) is assumed in the analysis at 
m n = 300 MeV, so high statistical calculations at more quark masses are necessary to obtain 
a definite conclusion for the quark mass dependence. We also leave this issue to studies in 
the future. 
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C. Comparison with other works 



In Fi g. [5] we compare our results (PACS-CS) obtained in 2 + 1 flavor QCD with those by 



ETMC [10l . and Lang et al. 12| in 2 flavor QCD. The upper panel shows the effective 



coupling constant g pnn and the lower panel displays the resonance mass m p as a function of 
m 2 . Here the systematic uncertainty for the determination of the lattice spacing is added 
to the statistical error in quadrature. A good agreement between our result and ETMC 
is observed for g p7r7r . The result for the coupling constant by Lang et al. takes a slightly 
smaller value, but it is almost consistent with other works. 

We see, however, large discrepancy for the resonance mass m p in the lower panel of Fig. |HJ 
One of possible reason for this discrepancy is the systematic error from the determination 
of the lattice spacing which is used to obtain m p and in the physical unit. In the present 
work the lattice spacing a = 0.0907(13) fm determined from ttiq in Ref. [l3| is used as 
explained before. ETMC used a = 0.079(2) (±2) fm given from the pion decay constant f n 



in Ref. [23]. In the work by Lang et a/., the authors determined it to be a = 0.1239(13) fm 
from the Sommer scale r = 0.48 fm as input. In order to avoid a spurious systematic error 
from the determination of lattice spacing, it is appropriate to compare our results with other 
works in terms of dimensionless quantities. In Fig. [9] we plot r m p as a function of (rom^) 2 
with the Sommer scale r . The value of r for the PACS-CS configurations has been reported 
as r Q /a = 5.427(51)(+81)(-2) [13] and that for ETMC as r /a = 5.32(5) H- In the figure 
the statistical error and the systematic uncertainty for the determination of r are added 
in quadrature. We see that the discrepancy between ours and ETMC tends to be smaller, 
but it still remains for the large quark mass. The result by Lang et al. takes a smaller 
value than those of the two works. The finite size effect can be considered as a possible 
reason of their small value of m p as commented by themselves in their paper. Their lattice 
extent La = 1.98 fm may not be large enough for their quark mass m n = 266 MeV. The 
three groups worked at a single lattice spacing, therefore the another possible reason of the 
discrepancy is the discretization error due to the finite lattice spacing. We can also consider 
several other reasons, the dynamical strange quark effect, the isospin breaking effect, the 
reliability of the parametrization of the scattering phase shift by fl30l and so on, but a 
definite conclusion can not be given here. A precise determination of the resonance mass m p 
by the calculation near or on the physical point closer to the continuum limit is an important 
work reserved for the future. 



IV. CONCLUSIONS 

We have reported on a calculation of the P-wave scattering phase shift for the isospin 
1 = 1 two-pion system and estimations of the resonance parameters of the p meson from 
the Nf = 2 + 1 full QCD configurations with a large lattice volume. The calculations are 
carried out at two quark masses, which correspond to = 410 MeV and 300 MeV. 

In order to extract the resonance parameters from the scattering phase shift, we 
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parametrize the resonant behavior of the P-wave phase shift in terms of the effective cou- 
pling constant g pn7T and the resonance mass m p . We find that this parametrization works 
well in the large energy region for our data at m n = 410 MeV and obtain g p7Tn = 5.52 ± 0.40. 

For m n = 300 MeV the statistics of our data is not enough to discuss the reliability of the 
parametrization. We leave an investigation on this point to the studies in the future. We 
carry out the fitting assuming that this parametrization also works in our energy region at 
m n = 300 MeV. Our result is g p7T7T = 5.98 ± 0.56, which agrees with the coupling constant 
at m w = 410 MeV and the experiment g p7T7T = 5.874 ± 0.014 within the statistical error. 
This suggests a weak quark mass dependence of the coupling constant. The studies at more 
quark masses are necessary to obtain a definite conclusion for the quark mass dependence, 
however. 

We find a discrepancy for the resonance mass m p among three lattice studies. Although 
a part of the discrepancy seems to be explained by different choices of the scale setting, 
other sources such as the discretization error due to the finite lattice spacing, the dynamical 
strange quark effect, the isospin breaking effect and the reliability of the parametrization of 
the scattering phase shift may be needed to resolve this discrepancy. Calculations near or 
on the physical point closer to the continuum limit are necessary for a precise determination 
of the resonance mass from lattice QCD. We leave this issue to studies in the future. 
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FIG. 1: ^fs/rrip for states tabulated in Table U on our gauge configurations at = 410 MeV 
(upper panel) and = 300 MeV (lower panel). 



0,ti 



G 



7V7V — >7V7V 



(*) 



-p,t 



V V 

0,t s p,t s 

p,t 






o, t s p,t s 




(Mi 



G 



p— >7T7T 




FIG. 2: Quark contractions of tttt — > tttt, tttt — > p and p —> irir components of the matrix of 
the time correlation function G(t). The time runs upward in the diagrams. The vertices refer to 
the pion or the p meson operator with the momentum at the time slice specified in the diagrams. 
The momentum p takes p = (27r/L)(0, 0, 1) for the A^" and p = (2tt/L)(1, 1, 0) for the 
representation. 
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t s = 12. The pion with the zero momentum in the sink operator (inr)(t) in (|10p is set at t\ = 42. 
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FIG. 8: Comparison of our results (PACS-CS) obtained in 2 + 1 flavor QCD with those by ETMC 
and Lang et al. in 2 flavor QCD. Upper panel shows the effective coupling constant g pn7T and lower 
is the resonance mass m p . The systematic uncertainty for the determination of the lattice spacing 
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Tables 



TABLE I: The ground and the first excited states with the isospin (I, I z ) = (1,0) for the irreducible 
representations considered in the present work, ignoring the hadron interactions. P is the total 
momentum, g is the rotational group in each momentum frame on the lattice and Y is the irreducible 
representation of the rotational group. The vectors in parentheses after (ir-ir) and p refer to the 
momenta of the two pions and the p meson in unit of 2ir/L. We use a notation (vr7r)(pi)(p2) = 
7r + (pi)7r - (p2) — 7r~(pi)7r + (p2) for the two-pion states. An index j for the representation takes 
j = 1,2,3 and k for the E~ takes k = 1,2. 



frame PL/(2n) g 
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( Pl + p 2 )(l,l,0) 



(vr7r)(e j )(-e j ) 

(7T7r)(e 3 + e fc )(-e fc ) - (vr7r)(e 3 - e fc )(e fc ) 
(7T7r)(0, 0,1) (0,0,0) 
(7T7r)(l,l,0)(0,0,0) 
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TABLE II: Results at m n = 410 MeV. P is the total momentum and T is the irreducible rep- 
resentation of the rotational group on the lattice. E is the energy extracted by fitting the time 
correlation function with the fitting range in a line of "Fit Range" . y/s is the invariant mass and 
k is the scattering momentum, which are related by y/s = y/E 2 — |P| 2 = 1\Jk 2 + m 2 . 5(k) is 
the P-wave scattering phase shift given by the finite size formulas in ([2]) . We use the value of the 
lattice spacing given in the previous work in Ref. [13J, a = 0.907(13) fm (1/a = 2.176(31) GeV), to 
obtain the values in the physical unit, where the error of the lattice spacing is not included. 
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/111 O/l > 7\ 

411.2(1. 7 ) 








frame 


CMF 


MF1 


MF1 


MF2 


PL/(2tt) 


(0,0,0) 


(0,0,1) 


(0,0,1) 


(1,1,0) 


r 


Tr 


E~ 


A 2 - 


Br 


Fit Range 


23-36 


23-36 


24-30 


24 - 30 24 - 30 24 - 30 


aE 


0.4200(49) 


0.4622(53) 


0.4426(24) 0.4774(27) 0.4891(30) 0.5364(60) 


CLyf! 


0.4200(49) 


0.4184(58) 


0.3967(27) 0.4352(30) 0.4027(37) 0.4589(70) 


(ak) 2 (xl0~ 3 ) 


8.4(1.1) 


8.1(1.2) 


3.63(41) 


11.63(57) 4.83(80) 16.9(1.5) 


y/s (MeV) 


914(11) 


911(13) 


863.2(5.9) 


946.9(6.4) 876.2(8.1) 999(15) 


s (GeV 2 ) 


0.835(19) 


0.829(23) 


0.745(10) 


0.897(12) 0.768(14) 0.997(30) 


a 2 k 3 / tan 5 /y/s (xl0~ 3 ) 


-1.872(19) 


-2.120(22) 5.79(85) 


-15.5(1.7) 2.36(49) -28(11) 


fc 3 /tan<5/Vi (xlO~ 2 GeV 2 ) 


-0.8865(88) 


-1.004(11) 2.74(40) 


-7.34(80) 1.11(23) -13.0(5.2) 



TABLE III: Same as Table H] for = 300 MeV. 



am n 0.1355(15) 
m,, (MeV) 294.9(3.3) 



frame 


CMF 


MF1 


MF1 




MF2 




PL/(2?r) 


(0,0,0) 


(0,0,1) 


(0,0,1) 




(1,1,0) 




r 


Tr 


E~ 


A 2 - 




Br 




Fitting Range 


23-36 


23-36 


24-30 


24-30 


24-30 


24-30 


aE 


0.400(13) 


0.452(12) 


0.3658(40) 


0.4511(85) 


0.4377(80) 


0.492(18) 


ay/s 


0.400(13) 


0.407(13) 


0.3086(47) 


0.4061(95) 


0.338(10) 


0.406(22) 


(ak) 2 (xlO" 3 ) 


21.7(2.5) 


23.0(2.7) 


5.45(59) 


22.9(1.9) 


10.3(1.9) 


22.8(4.5) 


^(MeV) 


871(28) 


885(29) 


672(10) 


884(21) 


736(23) 


883(48) 


s (GeV 2 ) 


0.758(48) 


0.783(50) 


0.451(14) 


0.781(37) 


0.542(33) 


0.780(84) 


a 2 k 3 /tan5/y/s (xl0~ 3 ) 


2.1(1.9) 


1.4(1.6) 


25(11) 


-6.60(97) 


4.1(4.3) 


-12.3(3.7) 


k 3 / tan 5 /y/s (xlO^GeV 2 ) 


1.00(90) 


0.68(76) 


12.0(5.4) 


-3.13(46) 


19(21) 


-5.8(1.7) 
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